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PREFACE 


This  technical  report  is  primarily  a  compilation  of  presentations  made  at  the  annual 
project  review  meeting  held  on  June  3-4, 1997  at  the  University  of  South  Carolina.  Additionally, 
this  report  includes  an  introduction  to  the  project,  and  copies  of  papers  that  were  written  under 
the  sponsorship  of  this  grant  during  the  past  year. 

Any  questions  regarding  the  content  of  this  report  or  specifics  of  the  project  can  be 
addressed  to 

Dr.  Roger  A.  Dougal 
VTB  Project  Director 

Dept  of  Electrical  and  Computer  Engineering 
University  of  South  Carolina 
Columbia,  SC  29208 

phone:  (803)777-7890 
fax:  (803)  777-8045 

email:  dougal@ece.sc.edu 
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IMTRODUCTIOfM 


The  Power  Electronic  Building  Block  (PEBB)  will  enable  the  Navy  to  meet  DOD  and  Navy  goals 
of  reduced  manning,  reduced  cost,  increased  effectiveness  and  enlianced  survivability  by  allowing 
radically  new  architectures  for  shipboard  power  systems.  Since  these  new  architectures  cannot  be  based 
on  historic  design  precedents  and  time-and-field-tested  design  rules,  extensive  prototyping  and  testing  are 
necessary.  These  prototypes  are  used  to  validate  the  designs  and  to  define  the  operational  envelope,  in 
both  intact  and  damaged  conditions.  The  use  of  virtual  prototypes,  rather  than  hardware  prototypes, 
allows  the  US  Navy  to  maintain  its  technological  superiority  by  exploiting  its  dominant  position  witli 
respect  to  information  technologies. 

The  Virtual  Test  Bed  (VTB)  provides  a  unique  capability  for  virtual  prototyping  of  PEBB  devices 
and  of  PEBB-based  electric  power  systems  by  integrating  into  a  single  simulation  environment  models 
that  have  been  produced  in  a  variety  of  simulation  languages  by  a  diversified  and  multi-technical  design 
team. 

Those  who  develop  new  pieces  of  the  shipboard  power  system  often  (and  for  good  reasons)  use 
different  software  tools  for  different  modeling  tasks.  Since  a  system  is  composed  of  many  such  entities, 
and  the  system  must  be  tested  as  a  whole,  one  encounters  the  need  to  compute  the  performance  of  a 
system  described  by  a  heterogeneous  collection  of  models.  The  primary  objective  of  the  VTB  project  is  to 
aeate  the  virtual  prototyping  environment  that  accepts  this  heterogeneous  collection  of  models  and 
integrates  them  into  a  single  simulation  so  that  a  design  engineer  can  evaluate  and  understand  the 
dynamic  performance  of  the  entire  system.  The  VTB 

•  allows  closer  collaboration  amongst  experts  in  different  fields 

•  allows  each  expert  to  use  the  best  design  tools  and  best  design  practices  within  their  own  area  of 

expertise 

•  allows  integration  of  more  aspects  of  the  design,  including  physical  configuration,  electrical 

configuration,  thermal  configuration,  etc.,  into  a  single  virtual  prototype 

•  eliminates  the  need  for  manual  translations  of  models  while  exploring  system  response 

•  allows  more  rapid  exploration  of  a  larger  design  space  to  yield  more  optimal  designs 

•  provides  more  advanced  visualizations  of  system  performance  to  help  build  an  intuitive 

understanding  of  the  influence  of  controllable  parameters 

The  approach  chosen  for  creating  the  VTB’s  mixed-language  virtual  prototyping  environment 
involves  translation  of  the  source  models  into  the  VTB  internal  language.  This  provides  a  flexible, 
powerful,  robust,  and  extensible  means  for  integrating  models  into  a  unified  simirlation  environment.  In 
addition  td  developing  the  model  translation  technologies,  the  VTB  project  also  advances  other 
technologies  associated  with  virtual  prototyping,  including  user  interfaces,  model  libraries,  execution 
environments,  and  visualization  tools.  These  technical  advances  allow  the  VTB  to  fully  support 
development  and  application  of  the  PEBB. 

A  secondary  objective  of  the  VTB  project  is  to  develop  a  base  library  of  models  that  can  be  used 
in  the  Virtual  Test  Bed.  These  models  should  execute  rapidly,  be  constructed  in  such  a  way  that  they  can 
be  connected  to  other  models,  and  yet  can  be  written  in  a  variety  of  modeling/simulation  languages. 
Particularly  important  to  the  study  of  PEBB-based  power  systems  is  an  accurate  model  of  a  PEBB. 

The  third  and  final  objective  of  this  project  is  to  support  development  and  applications  of  PEBBs 
by  using  existing  modeling  tools  and  existing  components  of  the  VTB  for  virtual  prototyping  of  the 
PEBB  before  the  VTB  is  fully  developed. 

The  presentation  materials  that  follow  describe  a  majority  of  the  technical  work  that  was 
performed  under  the  auspices  of  this  grant  during  the  past  year.  Time  constraints  during  the  meeting 
prevented  the  presentation  of  all  aspects  of  all  work  that  was  performed,  but  this  report  provides  a  fair 
cross  section  and  a  useful  summary. 

Somewhat  arbitrarily,  the  report  is  organized  according  to  the  presentation  scheme  that  was  used 
at  the  annual  review  meeting.  That  is,  items  are  arranged  according  to  the  site  at  which  work  was  done, 
rather  than  by  topic. 
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and  m  interface  nodes  to  the  Nonlinear 
Network.  Furthermore,  the  nonlinear  network 
has  n  internal  states. 


where  A3,  A4  are  nonlinear  functions. 
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k  is  the  Newton’s  iteration  count. 

E  is  the  mismatch  of  the  previous  iteration 
1  is  a  vector  with  ail  entries  equal  to  1 . 


Universal  Model  of  Components 

Procedural  Models  -  Example 


Turn-OFF:  Zero  Current  Crossing 


Universal  Model  of  Components 

Procedural  Models  -  Example 


If  in  ON-State  and  Zero  Current  Crossing  Occurred 


Network  Solver  -  Virtual  TestBed 
Interface 


Implementation  Issues 

Unified  Model  Representation 


(Matrix  Topology  vs  Network  Topology) 

Proven  Minimization  of  Multiply  Adds 


Dynamic  Range  (System  Stiffness) 
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Adaptive  Averaging  PEBB  Model  -  Independent  Time 
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Animation/Visualization  Engine 

Interface 
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Reliabilitv  Analysis 


Approximate  Analysis  of 
Permeable  Material 
Transmission  Lines 


E.  N.  Glytsis  and  A.  P.  Meliopoulos 


School  of  Electrical  and  Computer 
Engineering 

Georgia  institute  of  Technoiogy 
Atlanta,  Georgia 


OUTLINE 


Importance  of  Modeling  of 
Permeable  Materials 

Model  Description 

Example  Cases 

Summary  and  Future  Work 


Importance  of  Modeling 
Permeable  Materials 


Existence  of  Permeable  Materials  in  a 
Ship  Environment 

Nonlinear  Effects  on  Transmission 
Lines 

Magnetic  Field  Shielding 

Dependence  of  Series  Impedance  on 
Electric  Load 

Modeling  is  Necessary 

-  Safe  and  Stable  Electric  Power  System 

-  Accurate  Control  of  Electrical  Equipment 


Model  Description 


Transmission  Line  Equations 


oz 


m  r  T  ^[v] 
-r-  =  [c] 


oz 


ot 


Lumped  Element  Equivalent  Circuit 


[  V]  =  Az[r]  m  +  AzU]  4^  =  [R]  [/]  +  [L]  4^ 

ot 


ri  AT 

[j]  =  Az[c]-^  =  [C] 


ot 


ot 


ot 


Method  of  Conductor  Subdivision 


Subdivide  Conductors  Into  Segments 
Calculate  Geometric  Mean  Radii 
Calculate  Geometric  Mean  Distances 
Transmission  Line  Equations 

Vi  =l(l^j+Xy)ij 

j 

Combine  Segment  Equations  to 
Calculate  [R]  and  [L]  Matrices 

Self-Consistency  (iterative) 


Example  Case 

Copper  3/0  -  GRC  2” 
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Resistance  and  Reactance  Table 


Electric  Current 

(Amp) 

Rfor  100  ft 

(mOhms) 

X=aj:.  for  100ft 

(mOhms) 

600 

24.5232 

17.6632 

700 

23.2944 

16.8022 

800 

22.2744 

16.1369 

900 

21.3782 

15.5949 

1000 

20.5686 

15.1151 

1100 

19.8437 

14.6607 

1200 

19.2060 

14.2217 

1300 

18.6524 

13.7980 

1400 

18.1744 

13.3930 

Example  Case 

Copper  3/0  -  GRC  2” 


Magnetic  Field 


Relative  Permeability 


Example  Case 

Copper  3/0  -  GRC  2” 


Magnetic  Field  Density 


Summary 


Importance  of  Modeling  Transmission 
Lines  in  the  Presence  of  Permeable 
Materials 

Model  Description  Based  on 
Conductor  Subdivision  and  GMR, 
GMD 

Example  Cases 


Future  Work 


Calculation  of  Capacitance  Matrix 

Interface  Transmission  Line  Model  with 
Network  Solver 

Improvement  in  the  Calculation  of 
Magnetic  Field 

Comparison  of  the  Model  with  the 
Rigorous  Model  Based  on  Finite 
Element  Method 


Virtual  Test  bed  Project 

RIGOROUS  TRANSIENT  MODELING  OF 
PERMEABLE  MATERIAL  TRANSMISSION  LINES 


Karim  N.  Wassef 
Dr,  Andrew  F,  Peterson 

Georgia  Institute  of  Technology 

School  of  Electric  and  Computer  Engineering 


Introduction 


Why  pursue  this  study? 

The  need  to  simulate  signal  propagation  in 
transmission  lines  between  the  various  components  of 
the  Network  Solver. 

What  is  expected  of  the  model? 

The  model  should  be  capable  of  modeling  currents  and 
voltages  propagating  through  any  geometry  (expected 
to  consist  of  multiple  conductors  in  a  permeable  steel 
conduit) 

How  is  this  done? 

This  will  be  achieved  using  a  combination  of  the 
Finite  Element  and  Transmission  Line  Matrix  methods 


in  a  Time  Domain  formulation. 


Contents 


The  Transmission  Line  Matrix  Method 

Time  Domain  Formulation 

The  Need  for  Parameter  Matrices:  [R],  [L],  [C] 

The  Finite  Element  Method 

Magnetic  Vector  Potential  Formulation 
Obtaining  the  Parameter  Matrices:  [R],  [L],  [C] 
Graphical  User  Interface 

Computational  Requirements  and  Results 

3(|)  copper  lines  in  conductive  permeable  steel  conduit. 

Verification 

Analytical  ( Coaxial  Line ) 


Future  Work 


The  Transmission  Line  Matrix 

The  basic  equations: 

—V  =  -RT-L^T 

dz  dt 

^i=-3v-d^v 

dz  dt 

{ V }  Voltage  in  each  of  the  conductors 
{ I }  Current  in  each  of  the  conductors 

[R]  Line  Resistance 

[L]  Self  and  Mutual  Inductance 

[C]  Self  and  Mutual  Capacitance 

[G]  Line  Trans-conductance 

All  terms  are  per  unit  length  and  vary  as  a  function  of  both 
location  and  time.  The  formulation  uses  a  Finite 
Differencing  Scheme  (in  both  time  and  space). 


The  TLM  Matrices 


[L]:  per-unit-length  Inductance  Matrix: 


/i=/2=/3. 


=...=4=0 


[R]:  per-unit-length  Resistance  Matrix: 


[C]:  per-unit-length  Capacitance  Matrix: 


[C]  =  [P] 


-1 


G, 


Qx-Qi-Qi’  •  '—Qi-\-Qi+\-'  • 


The  Finite  Element  Method 


The  equation  to  be  solved: 


r 


V- 


-VA  \=J-<J 


dA 

dt 


M  j 

where  A  is  the  magnetic  vector  potential  and  Js  is  the 
forced  current  density.  The  material  properties  are: 


a:  conductivity 
|i:  permeability 

The  magnetic  field  vector  B  and  flux  <I>  can  be  determined 


as  functions  of  A: 


Ai  -4- 


A2  ◄- 


-H 


(g)  B 


^12 

/ 


—  A2  Aj 


Illustration:  Copper  Lines  in  a  Steel  Conduit 


Conduit:  ( Steel) 

Relative  Magnetic  permeability  |ir 
Conductivity  a  = 

Conductors:  (Copper) 

Relative  Magnetic  permeability  |ir 

Conductivity  c  = 

Excitation:  (Current) 

Number  of  Conductors:  3 

Amplitude:  1  Amp 

Frequency:  50  kHz 


1000 

8.3334E+06  S/m 

1 

5.9595E+07  S/m 


Phase: 


3  Phase 


Analytical  Verification 


L 

I 


In 


In 


ya  j 


+ 


h 


c^-b 


For  a=50  ,  b=100,  c=110  mm. 


Analytical:  L!  I  —  2.4846 X 10  H  !  t 

Finite  Element:  h!  i  —  2.4768x10  ^  H  !  I 


Future  Work 


Alternative  Vector  Formulation. 


Adaptive  Gridding  Algorithms. 


Model  Parameter  Interpolation. 


Methods  for  rapid  solution  of 


linear  equations. 


Graphical  User  Interface 


Mesh 


Magnetic  Vector  Potential  ( A ) 


Magnetic  Flux  Lines  ( A  contour ) 


Magnetic  Flux  Density  { I5I ) 


Vector  Magnetic  Flux  Density  ( B ) 


Computational  Requirements 

Approximate  number  of  unknowns: 

~  2000  nodes  ( 3500  elements  ) 
Matrix  density: 

~  0.004  (  highly  sparse  ) 

Number  of  iterations  /  time  step: 

n+1  (n  is  number  of  conductors) 
Number  of  spatial  iterations: 

i 

m  = - 


(I  is  length  of  the  transmission  line) 
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Future  Work  Planned 


Relation  of  Work  to  VTB  Pr 


-  High-Power  Digital  Switches 

-  Digital  Signal  Processor 


Relation  of  Work  to  VTB  Project 
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Exploits  symbolic  or  numerical  off-line  solution 
Leads  to  fast  and  accurate  on-line  simulation 


Relation  of  Work  to  VTB  Pr 
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Unknown  parameters 
Computational  efficiency 


New  Modeling  Results 
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New  Modeling  Results 


New  Modeling  Results 


New  Modeling  Results 


yk+i  =  C(dk)xk  +  D(dk)uk 


New  Modeling  Results 


VP  speiueo 


XPspeuueB 


New  Modeling  Results 


Time  (ms) 


New  Control  Design  Results 
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New  Control  Design  Results 
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-  Globally  stable  with  asymptotic  tracking 


New  Control  Design  Results 


New  Control  Design  Results 


New  Control  Design  Results 


New  Control  Design  Results 


Time 


New  Control  Design  Results 
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Time 


New  Control  Design  Results 


Time 


New  Control  Design  Results 
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nonlinear  repetitive  control 
wider  range  of  applications 


Future  Work  Planned 


Virtual  Test  Bed 


Two  Universities.  ‘Two  Aircraft  Companies. 

Biggest  in  the  Europe  Boiler  Factory.  *  Steel  mill. 

•  Several  electronic  companies. 


Virtual  Test  Bed 
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•  About  600  faculty  members  and  over  2500  staff 
members  and  research  assistants. 
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Virtual  Test  Bed 


fAi 


O 

^  ■— 
^  CL 

IP 


E  o 

O  (0 

?  o 
iS  o 
a. 

c  c 
o  o 

c  '"S  ■ 

■i  >  i2 

w  c 

i>-o  g 
0)0)  o 

.E  0  Q. 

“■§  £ 
ce  S  o 

o  ^ 

O  ■« 
Co® 
.=  ■“  0 

^  s  > 

0)-g  ® 

<D  ii 
^  O  C 

C  P  O 

■-  s-  0 
0  0  0 


o  n 


■5  i2  g 

^  H  (0 

w  <0  o 


(0 

75 

o 

3 

(0 

"O 

0 

> 

T3 

*o 

O 

c 

E 

■  HM 

0 

3 

TJ 

CQ 

H 

s 

ti 


vx 


b  5  -  ^  § 

^  o  'o  o  S  8 

^  <I  ^  ^ 

•  #  •  •  O  Tj 


Conclusion:  modular  approach  is  more  promising; 
if  both  can  be  used,  the  modular  way  is  preferable 
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Conclusion:  Unfortunately,  as  a  rule,  the  parts  in  power  system  are 
tightly  coupled.  The  use  of  modular  approach  may  cause  instability 
or  non-convergence. 

We  are  going  to  solve  this  problem. 


Virtual  Test  Bed 


Conclusion:  Unfortunately,  we  cannot  apply  these  methods  to 
improve  the  performance  of  modular  approach 
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We  propose:  multicycle  concurrent  and  We  propose:  modified  and  generalized 

iterative  approach  (next  reports);  coupling  approach; 


form  of  a  diagram  (see  next  slide). 


We  are  going  to  use  WR  (iterative  approach)  and 
CR  (one-sweep  approach); 
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Modified  coupling 


Virtual  Test  Bed 
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strong  global  feedback  as  well. 


Thermodynamical  models 
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types  of  excitations  from  environment  during  operation  (destabilizing 
factors)  that  allows  to  estimate  their  influence  on  the  circuit  operation. 
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The  movements  separation  method  in  modeling  of  power  systems 

Jordan  canonical  models 

A  processes  in  power  and  electrical  systems  can  be  described  in  form  of  matrix 
differential  equation  in  Couchy  form: 


Given  transfer  matrix  may  be  interpreted  as  model  of  parallel-cascade  structure, 
which  consists  from  separated  independent  submodels,  having  its  own  time 
constant  ^  >  where  -  i-s  eigenvalue  of  matrix  A. 

By  means  of  such  model,  electrical  processes  having  the  different  speeds 
should  be  located  in  different  part  of  the  considering  system. 

That  why  method,  based  on  transformation  of  initial  mathematical  model  to 
Jordan  canonical  form  we  call  as  movements  separation  method  or 
•  method  of  decomposition  in  the  time  domain. 


iiiwwiillHt'ilMf  ^ 


Virtual  Test  Bed 


The  structure  of  circuits  with  separated  movements 


1.  The  case  of  single-input  -  single-output  circuit: 


2.  The  case  of  multipoles  circuit 
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Partial  transfer  functions  which  corresponds 
to  all  possible  types  of  Jordan  cells 

From  structure  of  system  matrix  in  the  Jordan  basis  followed,  that  type  of  any  partial 
transfer  function  is  determine  by  type  and  multiplicity  of  correspondent  eigenvalue 
of  the  Jordan  cell  Hence,  for  obtaining  the  such  set  of  primitive  circuits  is  necessary  to 
investigate  all  cases  of  types  and  eigenvalues  multiplicity  of  system  matrix  A. 
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Simple  eigenvalues 


Eigenvalues  of  algebraic  multiplicity 


Eigenvalues  of  geometrical  multiplicity 


1.  5:i  m  Die  value: 


(p-A) 

3.  Simple  eigenvalue  of  algebraic 
multiplicity: 

K.(p)  =  X  I 


ni=l  k 


s  (p-^r 


5.  Simple  eigenvalue  of  geometric 
multiplicity: 


r. 


K.(p)  =  Ik.(p) 

j=l 


c'h' 

k^(p)  =  c;0;B:,  -  '  ■ 


p  - 


i.  Com.ra.eX”Coiiiiia;ate  pair  eigenvalues: 


a,p  +  a, 


^  ^  p^"l"bjP”i"bQ 

4.  Complex-conjugate  pair  eigenvalues  of 
algebraic  multiplicity: 

Jl  Ti^+l  „  n  4-  a 

‘‘OmkP  ^  "Imk 


K.(p)=SE 
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m  ’ 


6.  Complex-conjugate  pair  eigenvalues  of 
geometric  multiplicity: 

K,(p)  =  ik,(p) 

j=l 


ki(p) 


Hip  +  a, 


p^+b,p  +  bo 
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Realizations  of  all  possible  types  of  Jordan  cells: 


3-4.  Siriile  eigenvalue  or  complex-conjugate  pair  of  algebraic  multiplicityi  . 


5-6.  Simle  eigenvalue  or  complex-conjugate  pair  of  geometric  multiplicit/'  . 
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Reduction  of  circuits  with  separated  movements 


1,  Reduction  in  the  time  domain: 


In  the  Jordan  basis  system 
eigenvalues  are  separated. 
Due  to  any  eigenvalue 
determine  die  speed  of  partial 
transient  process  the 
movements  separation 
method  open  the  way  to 
system  reduction  in  time 
domain  by  removing  parts 
with  the  smallest 
eigenvalues. 


2.  Simultaneously  reduction  in  time  and  frequency  domains: 


I  1 2,73  -  2.80,  mod^l  re^uctiou  and  iraaisfoanadoTi 


In  the  Moore  basis  system 
singularvalues  are  separated. 
As  any  singularvalue 
determine  energetic 
contribution  to  the  system 
dynamic  behavior  then  we 
can  reduce  system 
simultaneously  in  time  and 
frequency  domains  by 
removing  parts  with  the 
smallest  singularvalues . 


By  means  of  (2.73-2.80) 
relations  between  Jordan  and 
Moore  basises  we  can 
accomplish  flexible 
reduction,  not  only  in  the 
time  domain,  but  also  in  time 
and  frequency  domains 
simultaneously 
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Controllability  and  observability  of  circuits  with  separated  movements 

. . . . 


CoiltrolkMeil  B,  AB, ... ,  aiB  linear  independent 


Observable  if  ...  C^-  ajg  bnear  independent. 


m  dse  Jordan  Basis 
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System  matrixincludes 
eigenvalues  of  algebraic 
multiplicity 


Controllability  and  obse¬ 
rvability  of  circuit  depends 
only  from  it  part  which 
doesn’t  include  eigenva¬ 
lues  of  algebraic  multipli¬ 
city. 


IP 


System  matrix  includes 

.eigenvalue  of  geometric 
multiplicity  ■ 

Circuit  is  uncontrollable  | 
and  unobservable  indepen¬ 
dently  of  controllability  and| 
observability  coefficients. 
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Conclusions 

•1 .  Transformation  of  operator  with  given  input-output  map  for  power  systems  and  its 
components  to  Jordan  form  is  equivalent  to  construction  of  system  with  such  structure 
that  processes  having  the  different  speeds  should  be  located  in  different  parts  of  that 
system.  In  this  case  we  have  system  with  time  and  space-separated  movements. 
Processes  in  different  parts  of  the  such  system  are  independent  one  from  another. 

•2.  Mathematician  model  of  any  electrical  network  in  Jordan  form  can  be  constructed  in 
form  of  parallel-sequential  structure.  In  the  general  case  this  structure  can  be  composed 
from  six  base  primitive  circuits  which  are  determining  by  type  and  multiphcity  of 
correspondent  eigenvalues  of  mathematician  model.  With  help  of  such  structure  the 
tasks  of  modeling  and  synthesis  can  be  solved  by  determination  number,  type  and 
parameters  of  primitive  circuits. 

•3.  Space  and  time  model  decomposition  gives  possibility  of  effective  using  analysis 
methods,  based  on  movements  separation  (direct  integration  methods,  adaptive  algorithms) 
as  well  as  algorithms,  based  on  space  decomposition  of  the  system  (relaxation  methods, 
iterative  algorithms  of  type  Gauss-Jacobi  and  Gauss-Zeidel). 

•4.  Jordan  canonical  models  are  very  effective  in  setting  the  hierarhial  families  of  circuit  and 
system  models.  They  are  effective  for  making  model  reduction  in  the  time  domain. 

•5.  For  simultaneous  reduction  in  the  time  and  frequency  domains  great  results  might  be 
received  by  transformation  the  Jordan  models  into  Moore  basis. 

•6.  The  general  features  (controllabiUty  and  observability)  of  circuits  with  such  models  have 
been  investigated. 

•V.The  circuit  synthesis  method  based  on  Jordan  canonical  models  have  been  developed. 

•8.  Suggested  model  is  an  invariant  for  a  set  of  power  (or  electrical)  systems  having  similar 
input-output  characteristics  and  can  be  used  as  abstract  class  for  models  derivation  in 
CAD  systems. 
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Electrical  circuits  synthesis  algorithm,  based  on  the  movements  separation 


Calculation  of  state  equations  for 
PEBB-ML  models 
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Pseudocode  of  the  algorithm; 

CStateEquations  StateEquations,  JordanEq; 

StateEquations  =  CalculateStateEquations(AFCArray  /*  TransientFunArray  */); 
JordanEq  =  TransformToJordan(StateEquations); 
if  (JordanBasisO) 

{ 

CListOfEigenValues  eig  =  JordanEq.GetEigenValuesO; 

CUstOfControl  ControlVector; 

CListOfObserve  ObserveVector; 

CJordanBlocksList  JordanBlocks  =  NULL; 

eig.FirstO; 

int  n_Multiplicity; 

for  (int  i  =  0;  i  <  eig.GetLengthO;  i-H-) 

{ 

if  (eig.GetTypeO  ==  TYPE_REAL) 

{ 

Multiplicity  =  eig[i].Ge(MultiplicityO; 
JordanBlocks.Add(GiveRealBlock(eig[i],  ControlVector[i], 

ObserveVector[i],  n_Multiplicity)); 

i  +=  n_Multiplicity  - 1; 

} 

else  //  Complex  pair  of  eigenvalues 

{ 

Multiplicity  =  eig[i].GetMultiplicityO; 
JordanBlocks.Add(GiveComplexBlock(eig[i],  ControlVector[i], 

ObserveVector[i],  n_Multiplicity)); 

i  +=  2*n_Multiplicity  - 1; 

} 

} 

ConnecJordanBlocks(&JordanBlocks); 

} 

else 

{ 

switch  (BASIS) 

{ 

case  PARALLEL: 

CalculateParallelRealization(JordanEq); 

break; 

case  CASCADE: 

CalculateSequentialRealization(JordanEq); 

break; 

case  MOORE: 

CalculateMooreRealization(JordanEq); 

break; 

case  FENCES: 

CalculateFencesRealization(JordanEq); 

break; 

default: 

break; 

} 

) 


wi  Jliiiffl  I 


.  1952  - 

^^lo  ENGINE^^ 


Virtual  Test  Bed 


An  example  of  circuit  realized  by  it  transfer  matrix 
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An  example  of  circuit  with  separated  movements 
which  equivalent  to  shown  on  previous  slide 
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Relations  between  the  movements  separation 
method  and  other  circuit  synthesis  methods. 


Simultaneously  reduction 
in  time  and  frequency 
domain  by  removing  die 
smallest  singularvalues 
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Chart  of  Transients  of  the  Conditions  Variables 
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5)  Future  Research 


COUPLED-CIRCUIT  MODEL 
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Number  of  branches  -  24 


PARK’S  CIRCUIT 


VOLTAGE-BEHIND-REACTANCE  MODEL 


Number  of  branches  - 10 


IN  LOAD 


Voltage-Behind-Reactance  Model  -12  s 


VOLTAGE-BEHIND-REACTANCE  MODEL 
(SATURATION  INCLUDED) 
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SABER  MODEL 
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Traditional  Detailed  Simulation 
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♦  Neglect  All  Circuit  Parasitics 

♦  Neglect  All  High  Order  Effects 


Need  For  High  Frequency 
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■  Mechanical  Failure  of  Bearings  Due  to 
Fluting 

■  Loss  Analysis 


H.F.  Simulation  Approach 


♦  Behavioral  Semiconductor  Models 

♦  Good  Parasitic  Models 


Uses  For  Proposed  Environment 
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■  Semiconductor  Device  Modeling 

♦  New  Behavioral  IGBT  Model 


Induction  Motor  Structure 
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Machine  Model 
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tandard  and  HF  Simulation  Mode 
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■  Actual  Response  Richer 
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♦  Behavioral  Models 


New  Behavioral  IGBT  Mode 


Conductance  Profiles 
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Recommendation 

—  Create  a  Visualization  Working  Group  from  the  VTB/PEBB  community  to  address  A  & 
V  issues. 
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Is  there  a  fundamental  reason  why  circuit-based  models  can  not  be 
coupled  to  differential  equation  based  models? 

—  It  is  an  issue  of  methodology.  Both  types  are  converted  to  some  kind  of  difference 
equation,  so  it  should  be  possible  to  connect  them. 


SOUTHCAROLINA 


iS 


C  0 
O 

s.i 

’5  0 

i 

o  § 

0  O 

-  _? 

"o  '*'* 

>>  0 
D)  0 
_o  C 

O  CO 
CL  C 
O  O 
Vj 

0  0 

JC  3 

O" 
C  0 

w  O 

.9  0 
t5  E 

■iZ  0 

-4— » 

0 

0  0 

®  f 

JD  0 

0  S 

^  M- 

0  o 

T3  ^ 

^  8. 

w  2 


■  ■PM 

c^ 

.o 

0 

P 

CO 

r 

;g 

P 

g 

•P 

To 

s 

1— H 

cd 

> 

a 

H- 

C/) 

2  o 

0  (0 
X)  Z3 
O 

CO  0 

E  ^ 

^  0 

0)  5 
x:  2 


'3  g 

•8  I 

'§  ^ 

«  o 


I  I' 

o 

00  (0 
^  > 
O  cd 


.jH  O 

'S  I* 
a 


c3  ^ 

^  I 

<0 

^  pH 


o  o 
PL,  00 


•p  c^* 
R  00 
ir^  <0 
O  Ph 

§  ^ 

^1: 

R  o 

pp  ^ 


(O  O* 

>  CO 

t>  (D 

^  g 

O  d 

rP 

00  CT* 


o  w 

’B  ^ 

CO  f  N 

a  -I 


9  ^ 


(D  O 

9  o 

a  « 

0  CO 

O  o 


*1H  (D 
CO 

^  p 
§  o 
K  P. 

r-H  ^ 

o  P 


SS  S' 

00  .a 


o 

o 

•I— » 

D) 

c 

'£ 

o 

o 

0 

-!-» 

13 

o 

JZ 


T3 

0 

0 

O 

Q. 

£ 

0 

JD 

C. 

0 

O 

0 

c 

_o 

V- » 

'l— 
>1— » 
0 
0 

>> 

c 

0 

£ 


o 

X 


c^- 

0 

£ 

o 

0 

c 

0 

•o 

3 


0 

0 

3 

O) 

c 

t: 

0 

D) 

>. 

X) 

0 

O) 

■D 

0 


O 

c 

c 

’0 

D) 

S' 

£ 

0 


O 

c 

0 

0 

c 

0 


3 

o 

ifc  *-■ 
“  0 

^  'S 
w  — 

•4— • 

0  ZJ 

ft 


CO 

o 

< 

0 

■D 

0 

■D 

C 

0 
■I— • 

0 

D) 

£ 

0 

x 

0 

0 

sz 

>4— » 

o 

XI 


0 

XI 

o 

D 

0 

0 


0 

> 

0 


O 

c 

o 

■D 

0 


W 

0 

■D 

O 

£ 

CQ 

h- 

> 


0 

■O 

0 

Ti 

C 

0 

-4-* 

0 

"0 


0 
Q. 

£ 
o 
o 

0 
T3 

0 

~o  O- 
£  ^ 
iS  0  "o 


c 

o 


0  0 
C  0 
0  ZJ 

D)  cr 

3  .2? 
O  £ 

0  O 


c 

'sz 


0 

X} 

0 

CO 

0' 

g 

'cl 

CO 


0 

0 

0 

C 

0 


0  O 

0  H= 

C  o 
0  0 

>  0 


SOUTHCAROLINA 


(0 

£ 


O  ‘tr 


0 

> 

0 


U) 

C 


0 

> 

0 


C 

0 


D 

Q. 

E 

o 

o 


0 

'+—  -  ■ 

o  C 

?  O 

?  £ 

o  2 

cO  . . 

•“  c 

U  0 

0  0 

Q  0 
00 

E 

_0 

JD 
O 

Q_ 


(Ji 

c 
o 
•<— » 

'0 

c 

0 


U) 

3 

O 


D) 

C 

'q. 

Q. 

0 
-I— • 

0 

I 

O 

0 

c 

0 


0 

3 

H-* 

0 

C 

C 


(0 

0 

E 

0 


-A— 

o  0 

0 

0  *“ 

^  E 

T5  0 


Will  collaborate  on  techniques  for  efficient/fast  computing 


Plans  for  next  year 


D 

•  rH 

> 

(D 

PQ 

H 

> 


June  1997 


SOUTHCAROLINA 


0)  0 


o  w 


D) 

0 


0 

CO  0 
O  ^ 

o  E 

Q.  ® 

O  £ 

0  C 

>  ^ 
m  0 

Q  ^ 


technologies  into  the  VTB  software 


SOUTHCAROLINA 


SOUTHCAROLINA 


SOUTHCAROLINA 


ll 


(fi 


o 

CO 

CO 

LLI 

o 

Q_ 

CO 

CO 

■  wmmm 

c 

LL 


Q. 

0 

O 

C 

o 

o 

0 

O 

0 


O 

0 

0 

c 

0 


0 

iD 

0 

CO 


D) 

0 

00 


0 

0 

U 

O 

£ 

0 


O 

£ 

0 


■D 

o 

TO 

■  HHI 

0 

£ 

o 

c 

o 

■ 

o 

3 

0 

O 

0 

0) 

c 

T3 

0 

0 

TO 

T3 

0 

0 

U 

0 

o 

O 

■ 

c 

o 

0 

mmJ 

O 

£ 

c 

o 

0 

Q. 

O 

£ 

0 

c 

■  ■■Hi 

o 

> 

TO 

0 

0 

0 

Q 

CD 

Q. 

0 

0 


0 

0 

CD 


o 

c 


O 

0 


0 

0 

0 

0 

0 

0 

TO 

0 


VNnoavDHxaos 


CO 

D) 

b) 

c 

o 

■ 

■D 

Q. 

13 

c 

O 

■  ai^Hi 

C 

CO 

■  HHBHi 

E 

o 

CO 

c 

■ 

■ 

b 

> 

C 

0 

o 

C 

CO 

o 

0 

Q. 

E 

o 

b) 

o 

0 

Q. 

c 

■  MMi 

O 

0 

3 

0 

LL 

Q 

• 

• 

CO 


ID 


0 

O) 

0 

C 

0 

E 

0 

0 

0 

-Q 

0 

0 

TD 


^  0 


0 

C 

o 

Q. 

E 

o 

o 

O) 

■ 

;g 

■ 

13 

I 


0 


O) 


O 

0 

> 

0 

TD 

I 


SOUTHCAROLINA 


U-! 


CD 

CD 

LU 

CL 

'o 

0 

■D 

O 

E 

T3 

O 

O 

D) 


0 

O 


O 

O 

D) 


O 

0 

0  O 

Q  .E 


"D 

O 

0 

£ 

o 

0 

0 

E 

> 

to 

0 

0 

> 

c 

0 

■  mamma 

o 

■  ammam 

0 

Q. 

0 

0 

Q. 

O) 

O 

C 

0 

■  MM 

c 

0 

0 

Q 

0 

0 

O 
■  ■■■■ 

> 

0 

TJ 

■D 

0 

O 

C 

0 

> 

■D 

0 

O 

0 


0 

E 

Q. 

O 

0 

> 

0 


O) 

'(/} 

0 

■D 

CD 

CD 

LU 

CL 


0 


O 

E  _  E 
clO  ^ 

-2  0)  J 
0  .  0 

Q  >2.  < 


